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Many body localization in the presence of a single particle mobility edge 
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In one dimension, noninteracting particles can undergo a localization-delocalization transition 
in a quasiperiodic potential. Recent studies have suggested that this transition transforms into a 
Many-Body Localization (MBL) transition upon the introduction of interactions. It has also been 
shown that mobility edges can appear in the single particle spectrum for certain types of quasiperi¬ 
odic potentials. Here we investigate the effect of interactions in two models with such mobility 
edges. Employing the technique of exact diagonalization for finite-sized systems, we calculate the 
level spacing distribution, time evolution of entanglement entropy, optical conductivity and return 
probability to detect MBL. We find that MBL does indeed occur in one of the two models we study 
but the entanglement appears to grow faster than logarithmically with time unlike in other MBL 
systems. 
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Introduction: Noninteracting particles in the pres¬ 
ence of disorder exhibit the phenomenon of Anderson 
localization [l]. In one and two dimensions an arbi¬ 
trarily weak amount of disorder is sufficient to localize 
all eigenstates 3- I n three dimension a mobility 
edge, defined as a threshold eigenstate with energy E c 
that separates localized and delocalized states can ex¬ 
ist. The question of how Anderson localization is mod¬ 
ified in the presence of interactions has become an area 
of intense activity following the seminal work of Basko, 
Aleiner, and Altshuler [dj. These authors argued that an 
interacting many-body system can undergo a so called 
Many-Body Localization(MBL) transition in the pres¬ 
ence of quenched disorder. This MBL transition involves 
highly excited many-body quantum states and can thus 
extend up to even infinite temperature in contrast to a 
usual quantum phase transition [E], which involves only 
the ground state. Traditional notions of statistical me¬ 
chanics do not apply to this transition and the localized 
phase, including the Eigenstate Thermalization Hypoth¬ 
esis (ETH) [6H8[ for the mechanism of thermalization in 
isolated quantum systems SO- It has thus been sug¬ 
gested that there are emergent conservation laws for these 
localized systems 11 , |12i[ like for integrable ones, which 
too do not thermalize [ 13 , 14]. 

It is possible to have a localization-delocalization tran¬ 
sition similar to the MBL transition for non-interacting 
one dimensional models with quasi-periodic potentials. 
An example of such a model is the Aubry-Andre 
model [H|(AA model) , which has the form 


H = ^2 hiUi - t (cjcj+i + ct +1 Cj) , s 


(1) 


where c (c^) annihilates (creates) spinless fermions, t is 
the hopping and hi an onsite potential with the quasi¬ 
periodic form hi = /icos(2-7rm + <f >), where a is an ir¬ 
rational number and <f> an offset. This model has a 
localization-delocalization transition at h = 2t, where 
all states are (de)localized for h(<) > 2 1. A numerical 


study of this model with a nearest neighbor interaction of 
the form V V ; rqrq+i has shown that the single particle 
transition changes into an MBL transition akin to the 
one in models with on-site disorder [l6j]. Furthermore, 
this model has recently been emulated in experiments on 
cold-atoms in the non-interac ting limit [171 Il8j and with 
interactions to observe MBL [19[ . 

Modifications to the AA model have been proposed 
to yield models which possess single-particle mobility 
edges [ 20 M 22 I ]. It has been argued that in the presence of 
(even weak) interactions, localized single particle states 
can thermalize when coupled to a bath of even a single 
delocalized state that is protected topologically or oth¬ 
erwise [231 ]. However, a very recent work shows that in 


the absence of any such protection, under certain condi¬ 
tions, the localized states can instead thermalize the bath 


24]. Models with single particle mobility edges are ideal 


to study the latter scenario since the delocalized states 
have no protection from localization. 


In this work, we study whether MBL occurs in mod¬ 
els with single particle mobility edges upon switching on 
weak to moderately strong interactions. Employing ex¬ 
act diagonalization of finite-sized systems, we calculate 
various diagnostics to detect MBL such as the level spac¬ 
ing distribution, time evolution of entanglement entropy, 
optical conductivity and return probability. Our conclu¬ 
sion is that MBL occurs in one of the models we study 
but not the other. However in the localized phase we 
observe, the entanglement entropy increases appears to 
increase linearly with time (like in an ergodic phase) but 
saturates to a sub-thermal value characteristic of MBL. 
The growth of entanglement entropy with time in a regu¬ 
lar many-body localized phase is logarithmic [iEl |2 t 1 . All 
other diagnostics appear to be consistent with regular 
MBL. We examine the possible reasons for the different 
behaviors of the two models and also provide possible 
reasons for the observed linear growth of entanglement 
with time. 


We have studied two different interacting one- 
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dimensional models of spinless fermions, which in the 
non-interacting limit have single particle mobility edges. 
The first, which we shall refer to as model I is described 
by the Hamiltonian 

H = ^2 hiTii - t(c\c i+ 1 + cj +1 Cj) + Vriin i+1 (2) 

i 

where hi = /icos(27rm n + <j)) with 0 < n < 1. For V = 0 
and n= 1, this is just the AA model. However, for n < 1 
and V = 0, the model has a single-particle mobility edge 
when h < 2t [io|, (2lj. All single particle states with 
energy between ±|2t — h\ are delocalized and all other 
states are localized. For h > 2t all single particle states 
are localized as in the usual AA model. 

The other model (which we refer to as model II) is also 
of the form in Eqn. [2] but with 

1 — cos(27T7o: + < f >) 

1 1 + /3 cos(2nia + <j>) 1 

with (3 £ (—1,1). When (3 = 0 and V = 0, this model 
also reduces to the AA model. For V = 0, there is a 
mobility edge separating, localized and extended states 
at an energy E given by f3E = 2 (t — h/2). This model 
also can be experimentally realized [22]. 

We have studied both models using exact diagonaliza- 
tion on finite-sized systems up to size L = 16 (data in 
plots shown for L = 14) with open boundaries and have 
averaged over the offset </> for better statistics, t = 1 and 
a = i n all our calculations. We now discuss our 

results. 



FIG. 1: (Color Online)The variation of the mean of the 
ratio between adjacent gaps in the spectrum for L = 14 
at half filling for model I and model II. The blue dotted 
line is for the Poissonian distribution and the pink one is 
for the Wigner-Dyson distribution. 

Energy level spacing statistics: Energy level spacing 
statistics is often used to characterize the MBL tran¬ 
sition. There is a crossover from a Wigner-Dyson to 
Poissonian distribution upon going from the ergodic to 
many-body localized phase, which can be tracked by the 
ratio of successive gaps, r n = where 



FIG. 2: (Color Online)Variation of the Renyi entropy 
for L = 14 at half filling for the two models with 
different parameters. 


S n = E n+ 1 — E n , the difference in energy between the 
n th and n + 1 st energy eigenvalues. For a Poissonian 
(Wigner-Dyson, specifically of the Gaussian Orthogonal 
type) distribution, the mean value of r is 2 In 2—1 ss 0.386 
(« 0.5295). The distribution function P(r) —7 0, as 
r — > 0 in the presence of level repulsion. 

For model I, with V = 0 , h < 2 and n = 1, all single 
particle states are delocalized. As V is increased, the 
level spacing distribution starts to follow the Wigner- 
Dyson distribution. For, n < 1, with a mobility edge, 
level statistics obey the Wigner-Dyson distribution, even 
though there are localized states as shown in Fig. [l] Deep 
in the localized phase (h » 2), increasing V yields a 
Poissonian distribution in both cases (n = 1.0 and n < 
1 . 0 ) 

Unlike for model I, the position of the mobility edge 
in the non-interacting limit of model II can be tuned 
by varying the parameters /3 and h We choose, 

h = 8 and change [3 from -0.95 to 0 so the fraction of 
single particle localized states increases progressively. In 
contrast to model I, here the level spacing distribution 
appears to be Poissonian for V ^ 0 as can be seen in 

Fig. CD 

Entanglement entropy: The entanglement entropy is 
another diagnostic that can be used to distinguish be¬ 
tween the ergodic and many-body localized phases. We 
have studied the time evolution of the entropy S(t) by 
sampling the initial unentangled states at random over 
the entire energy spectrum, which is equivalent to work¬ 
ing at infinite temperature [l^ ■ S(t) has been argued to 
grow linearly in the ergodic phase and logarithmically in 
the many-body localized phase (25], [28| . 

The system of length L is divided into two equal parts 
A and B. Our calculation is of the order 2 Renyi en¬ 
tropy i5>2(t) = — log 2 {TrAPA{t) 2 ) (which is computation¬ 
ally less expensive than the von-Neumann entropy) j29| , 
where pA{t) is the reduced density matrix of A obtained 
from the instantaneous state of the full system. It is 
known that in the ergodic phase, S 2 (t) ~ i at long times 
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FIG. 3: (Color Online)Thc variation of 
SS = S 2 (t, V ) — S 2 (t, V = 0). S 2 is the Renyi entropy 
for L = 14 at half filling for the two models with 
different parameters. The dotted lines are linear fits in 
t. (Inset) The variation of the saturation value of S 2 
with L. The blue dotted line corresponds to thermal 
value of *5*2 — — 1.2 for system size L. 


and saturates to the infinite temperature thermal value 
while for the usual many-body localized phase with weak 
interactions, S 2 (t) ~ £log(i), where £ is the localization 
length of the single particle eigenstates. It saturates to a 
value much smaller than the thermal value, but which is 
still extensive in system size. For our system, the infinite 
temperature 5*2 ~ 4 — 1.2 for system size L [16|. 

For model I with a single-particle mobility edge, S 2 (t) 
increases linearly with time but then appears to saturate 
to the thermal value as shown in Fig[2] However, for 
model II, S 2 appears to grow linearly with time but sat¬ 
urates to a value smaller than the thermal value. This 
can be seen from Fig. (2j where the time evolution of S 2 
has been plotted for model II for V = 0.2 , h = 8 and 
/3 = —0.95 ,—0.75 and —0.6, with progressively increasing 
fractions of single-particle localized states. The satura¬ 
tion value depends on the number of localized single par¬ 
ticle states: As the fraction of single-particle delocalized 
states increases, so does the saturation value. 

To confirm the linear growth S(t), we have plotted 
SS = S 2 (t,V ) — S 2 (t,V = 0) in Fig. [3] as a function of 
time. At very early times S 2 (t, V ) and S 2 (t, V = 0) tend 
to coincide, reflecting the formation of short range en¬ 
tanglement at the cut between the subsystems. Then, 
S 2 (t, V = 0) saturates but for the interacting system, S 2 
keeps growing with time as shown in Fig. [3] At intermedi¬ 
ate times, as long as there is a mobility edge in the single 
particle spectrum, SS fits quite well to a linear function 
of time. When all single particle states are localized, the 
growth of SS as a function of t is much slower than linear 
and possibly logarithmic. At long times, SS saturates to 
a sub-thermal value in all cases (For a calculation to even 




FIG. 4: (Color Online) (A)The variation of a(u) with w 
for the two models for L = 14 at half filling for V = 0.8. 
The rescaled values of a are plotted . 10cr, er and O.ltr 
are plotted respectively for model I (n = 1, ft = 1.5 and 
n = 0.4, h = 1.5) and model II (h = 8, /? = —0.8). The 
dashed lines are the best fit lines. (B)The variation of 
the return probability C(t ) as a function t for different 
models for L = 14 at half filling. The dashed lines are 
the best fit lines 


longer times see the supplementary information HI). 

We have also plotted the saturation value of S 2 as a 
function of system size L. As shown in the inset of Fig. OH 
~ L for the ergodic phase as well as for the model 
with a mobility edge. This plot also shows that the S ^ 0,1 
curve for the system with the single particle mobility edge 
system does not intersect the curve for the ergodic system 
when extrapolated to the thermodynamic limit. Thus, 
the saturation of the entropy to a sub-thermal value is 
not a finite-size effect. 

Optical conductivity: The optical conductivity cr(u>) 
is another diagnostic that can be used to identify the 
ergodic and many-body localized phases. In the case of 
a clean metal, the DC conductivity a(u> = 0) ^ 0 with 
a frequency-dependent additive term that goes as w 1 / 2 
at high temperature (&!}. In the presence of disorder, a 
subdiffusive phase can exist even on the thermal side of 
the MBL transition [32l435| . for which cr(u>) ~ to a with 
0 < a < 1. In the many-body localized phase, cr(u) ~ w“ 
with 1 < a < 2 [36J and a —> 1 as the transition is 

approached, cr(uj) is given by the Kubo formula, 

Ta(ui) = I < > 1 2 S(uj-E m + E n ) (3) 

ran i 

as T —>• oo, where, m,n are the many body eigenstates 
of the system with energies E m and E n . ji is the local 
current density. 

As shown in Fig. SKA) for model I at very low values 
of w, cr(w) ~ uj 3 . This is from a combination of level 
repulsion (one power of w) and open boundary conditions 
(w 2 ) Q. Subtracting this out, we obtain a ~ w 1 / 2 and 
a ~ a ; 3 / 4 for model I with h = 1.5, n = 1.0 and h = 1.5 
and n = 0.4 respectively. For model II in the presence 
of a single-particle mobility edge, after subtracting out 
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the uj 2 dependence [37J, a(oj) ~ w“ at low frequencies 
with 1 < a < 2 like in the usual many-body localized 
phase. In Fig. QJA) for a particular choice of parameter 
(3, a ~ cu 3 / 2 . 

We have verified that we obtain the same exponent a 
even with periodic boundary conditions, where the sub¬ 
traction is of a different power of oj. Further, we find 
that the exponent for model II increases as the fraction 
of localized states for V = 0 increases consistent with the 
expectation that the system gets pushed deeper into the 
many-body localized phase if it starts with more localized 
states without interactions. 

Return probability: The return probability C(t), mea¬ 
sures the probability of particles to return to their initial 
positions during the evolution of the system and is de¬ 
fined as 




, 'i'bJmn 1 1 


< n|(rij — 1/2 )|to > 


( 4 ) 


where, Z is the Hilbert space dimension. We have cal¬ 
culated C(t) = j^^2jC3{t) with C(t = 0) = 1. In the 
ergodic (diffusive) phase, C(t) ~ t -1 / 2 and in the many- 
body localized phase, it remains finite in the long time 
limit The behaviour of C(t) at long times is dras¬ 

tically different for the two models as can be seen in 
Fig. SKB) . For model I, with n = 0.4 and h = 1.5 at 
long times, C(t) ~ t~ b with b = 0.18 and for model II 
with /3 = —0.95 and h = 8, C(t) does not decay with 
time. The result for model I is consistent with the scal¬ 


ing relation a + 2b = 1, proposed by Agarwal et al. [34 [. 
We note however that we have not been able to clearly 
observe C(t) ~ f -1 / 2 in the thermal phase probably due 
to limitations of system size. 

Discussion: We have demonstrated the effect of in¬ 
teractions on models with mobility edges in the non¬ 
interacting limit. Our numerical results employing a 
number of different diagnostics show that an MBL phase 
can occur in such a situation. We find that model II dis¬ 
plays MBL while model I does not. One possible reason 
for this is finite size effects. To examine this, we have 
calculated the Inverse Participation Ratio (IPR) for all 
the states of an L = 14 system with V = 0 

The IPR of a normalized eigenstate i/j is defined 
IPRq, = )T/ ; | cj | 4 , where Cj is the amplitude of i/j at site 
j. IPR ~ 1 for a localized state and is much smaller 
(typical ~ 1 /L) for a delocalized one. The IPR val¬ 
ues for the two models are shown in Fig. [5] It can be 
seen that while there are localized states (with IPR of 
order 1) along with delocalized ones for model II, the 
states of model I appear to be delocalized for our system 
size. This behavior presumably persists even with in¬ 
teractions (which generally tend to cause delocalization). 
As a result, none of the diagnostics for this model show 
any evidence of many-body localization, to observe which 
probably requires larger system sizes that are not easily 



FIG. 5: (Color Online) IPR as a function of energy 
levels for the two models for L = 14. 


accessible with exact diagonalization. Another possibil¬ 
ity is that even in the thermodynamic limit the localized 
states of Model I are only “weakly” localized compared 
to those of Model II and thus fail to localize the bath. 
Another possibility is that the delocalized states of model 
I are “inherently” more robust compared to the localized 
ones whereas for model II, it is the other way around. 
A calculation of the IPR for large system sizes seems to 
suggest that this is true 0]. Additional calculations of 
the matrix elements for flip-flop processes and the lo¬ 
calization length exponent v[? } seem to indicate that 
model II may not satisfy the conditions for delocaliza¬ 
tion even when coupled to a protected bath [23J. Thus, 
an introduction of interactions would tend to cause MBL 
in model II and thermalization in model I even in the 
thermodynamic limit. Additional studies are required to 
fully understand the differences between the two models. 

For model II, the entanglement entropy appears to 
grow linearly with time (instead of logarithmically) be¬ 
fore saturating to a sub-thermal value. A possible ex¬ 
planation is the simultaneous but independent contribu¬ 
tions of the delocalized and localized states which indi¬ 
vidually would produce linear and logarithmic growth re¬ 
spectively. For sufficiently long times, the linear growth 
would dominate, which is what we observe. A mech¬ 
anism has been proposed recently invoking the idea of 
rare thermal regions in a many-body localized phase jjjjijj 
mainly to explain the behavior of <j{uf) near the MBL 
transition. The specific systems studied has spatially sep¬ 
arated ergodic and thermal regions. A calculation of S(t) 
performed by us for a similar system yields faster than 
logarithmic (algebraic) growth of S 2 with time, similar 
to what we observe for model II. Thus, the delocalized 
states in our models could be performing a role anal¬ 
ogous to that of rare thermal regions and producing a 
linear growth of entanglement. Note however, that the 
quasi-periodic potential in our models is correlated at dif¬ 
ferent sites and so no true rare regions in the sense of [sfij 
can actually occur. The analogy is therefore not a deep 
one. Algebraic growth of entanglement with time in a 
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many-body localized phase has also been observed in the 
presence of long-range interactions [381 ]. Another possi¬ 


bility is that the apparent linear growth is a finite-size 
effect in model II and will eventually become logarithmic 
for sufficiently large system sizes. The slow growth of the 
localization length near the mobility edge for model II as 
characterized by the exponent v [30| might be a sign that 
finite-size effects are important. 

Note added: A related study of many-body localiza¬ 
tion in model II appeared at the same time as ours [iffli ] 
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